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We study critical properties of the relaxation time at a threshold point in switching processes be- 
tween bistable states under change of external fields. In particular, we investigate the relaxation pro- 
cesses near the spinodal point of the infinitely long-range interaction model (the Husimi-Temperley 
model) by analyzing the scaling properties of the corresponding Fokker-Planck equation. We also 
confirm the obtained scaling relations by direct numerical solution of the original master equation, 
and by kinetic Monte Carlo simulation of the stochastic decay process. In particular, we study the 
asymptotic forms of the divergence of the relaxation time near the spinodal point, and reexamine 
its scaling properties. We further extend the analysis to transient critical phenomena such as a 
threshold behavior with diverging switching time under a general external driving perturbation. 
This models photo-excitation processes in spin-crossover materials. In the ongoing development of 
nano-size fabrication, such size-dependence of switching processes should be an important issue, and 
the properties obtained here will be applicable to a wide range of physical processes. 
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I. INTRODUCTION 

Relaxation phenomena in strongly interacting systems 
have been studied extensively, including various types of 
threshold phenomena. One of the most typical examples 
is the decomposition at the coercive field (the end of the 
hysteresis loop of a ferromagnetic system). This phe- 
nomenon appears in the field dependence of the order 
parameter in the ordered phase. In mean-field studies, 
this phenomenon is well described by the change in the 
free energy as a function of the magnetization (the order 
parameter of a ferromagnetic system). That is, in the or- 
dered state, the free energy has two minima representing 
the symmetry-broken states. When we apply an external 
field, one of them is selected to be the equilibrium state. 
For a weak field, however, the other state remains as a 
metastable state. When the field becomes strong, the 
metastable state finally becomes unstable. This point 
is called the spinodal point. However, in models with 
short-range interactions, fluctuations cause the system 
to escape from the metastable state through nucleation 
phenomena, and thus the change of the relaxation time is 
smeared. Thus, although there are several studies on the 
divergence of the relaxation near the spinodal points [1] , 
the spinodal phenomenon in short-range systems must 
be defined only as a crossover [2] . 

Recently, however, it has been pointed out that the 
mean-field universality class is realized in spin-crossover 
type systems [3]. In these materials, the volume of the 
unit cell changes, depending on the local two-fold states 
(say the high-spin (HS) and low-spin (LS) states). The 
volume change causes a lattice distortion, and the elastic 
interactions among the local distortions cause an effective 
long-range interaction among the spin states. The criti- 



cal properties of the spin-crossover system turned out to 
belong to the mean-field universality class. It was also 
found that the finite-size analysis for various quantities is 
very similar to the long-range, weakly interacting model 
(the Husimi-Temperley model). 

In spin-crossover systems, one may expect that the dy- 
namics also corresponds to that of the mean-field model. 
In particular, various threshold phenomena have been 
pointed out in the dynamics of spin-crossover type mate- 
rials. For example, a spinodal phenomenon without nu- 
cleation type clustering was reported in numerical calcu- 
lations [3] . The change is considered not to be a crossover 
process but a change with a true critical singularity that 
can be described by mean-field dynamics [4]. Moreover, 
threshold phenomena have been pointed out in the de- 
pendence of the photo-excitation process on the strength 
of the photo-irradiation, which can be modeled as a kind 
of spinodal phenomenon [5]. The metastable state does 
not relax to the stable state in the mean-field approxi- 
mation, which corresponds to the infinite system size in 
the Husimi-Temperley model. 

In recent extensive studies on nano-size systems, finite- 
size effects turn out to play important roles. Therefore, 
it would be useful to study finite-size effects on the spin- 
odal phenomenon as a critical dynamical process. The 
critical phenomena of the phase transition were studied 
in a dynamical mean-field model. The relaxation time 
diverges when the parameter approaches the threshold 
value. Binder studied the divergence of the relaxation 
time in spinodal phenomena in the mean-filed model by 
a Monte Carlo method [1]. Although the metastable state 
does not relax to the stable state in the mean-field ap- 
proximation, Paul, Heermann, and Binder studied the 
relaxation from the metastable side in finite-size systems, 
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and obtained a finite-size scaling form of the relaxation 
time [6] . The size dependence of the relaxation was also 
discussed recently [7] . 

In the present paper, we study the asymptotic behav- 
ior of the relaxation time of the infinitely long-range in- 
teraction model (the Husimi-Temperley model) near the 
spinodal point including the metastable region, the spin- 
odal point, and the unstable region. We reexamine the 
scaling properties near the spinodal point, which have 
been proposed by Paul, Heermann, and Binder [6]. 

As in the previous studies, we adopt the Glauber dy- 
namics, and derive a master equation for the probability 
distribution of the total magnetization. We first derive a 
master equation as a function of the total magnetization, 
which is possible in the Husimi-Temperley model because 
of the long-range nature of the interactions. Then, we de- 
rive a Fokker-Planck equation by using an expansion in 
terms of the inverse system size, which is an example 
of the van Kampcn Omega-expansion [8, 9]. We ana- 
lyze the Fokker-Planck equation and derive a scaling re- 
lation and also asymptotic forms of the relaxation times. 
These properties are confirmed by direct numerical in- 
vestigations of the original master equation, as well as 
corresponding Monte Carlo simulations of the stochastic 
decay process. 

We also extend the analysis to general threshold phe- 
nomena, such as switching from LS to HS states in spin- 
crossover materials by photo-irradiation. We obtain the 
effects of the photo-irradiation on the master equation, 
and show that the critical properties of these processes 
are the same as those of the spinodal phenomena. 

The rest of this paper is organized as follows. In Sec- 
tion II, we review the spinodal phenomena and define 
the master equation for the long-range model. Then, we 
study the asymptotic size dependence of the relaxation 
time near the spinodal point in Section III. Numerical 
confirmation of the asymptotic forms is given in Section 
IV. Next, we extend our study to threshold phenomena 
under external pumping in Section V. Finally, a summary 
and discussion are given in Section VI. 



II. INFINITE LONG-RANGE MODEL AND 
THE SPINODAL POINT 

We investigate the relaxation phenomena near the 
spinodal point in the infinitely long-range model. This 
is a spin model in which each spin interacts equally with 
every other spin. The Hamiltonian is 



H = 



_J_ 

'2N 



N 



M 



HM, M = J2 a i 



(1) 



where H is the magnetic field and <7j = ±1. It is well 
known that the mean-field theory is exact for this model 
in the limit of N — > oo. This model shows a second- 
order phase transition at T — Tq = J, H = 0. Below the 
critical temperature, a metastable state exists for weak 



magnetic fields. When the magnetic field becomes strong, 
the metastable state becomes unstable at a certain point, 
known as the spinodal point. In order to determine the 
spinodal point, we consider the extended free energy, i.e., 
the free energy for fixed total magnetization. In the in- 
finitely long-range model, this is given by 

/(to) = -^m 2 - Hm - In A rC( A r +M )/ 2 



1 / 1 + TO 1 + TO 1 — TO 1— TO 



(2) 



in the limit N — > oo where we use Stirling's formula 
and m is the magnetization per spin (to = M/N). The 
spinodal point is given by the following conditions, 



5/ 
dm 



and 



dm 2 



0, 



which give 



1 i ' V fl7 

at which the magnetization is given by 



m S p = ±\ 1 - 



(3) 



(4) 



In this paper, we consider the case where we increase the 
field from a negative value. Therefore, we consider the 
behavior of a locally stable state at negative magnetiza- 
tion around tosp < 0. 

We study the dynamics via the standard master equa- 
tion 



dP(S, t) 
dt 



J2Ws^s'P(S) + J2 w s'^sP(S r ), (5) 



where S and 5" denote states of the system and Ws^s' 
is a transition rate from S to S'. The probability of the 
state S at time t is denoted by P(S,t). 

Among the many possible dynamical models (choices 
of the transition rate) , we adopt the Glauber dynamics in 
this work. In the Glauber dynamics, the transition takes 
place as a flip of a local spin, and the transition rate Wij 
from a local spin state i to a local spin state j is given 
by 



Wij = 



T e~> 3Ei + e~P E J 



(6) 



where Ei denotes the energy of the system in spin state 
i, and To is some characteristic time scale. In this paper, 
we scale the time by To and set To = 1 for simplicity. 

With this transition rate, we construct a master equa- 
tion for the mean-field model. As the Hamiltonian de- 
pends only on the magnetization M, the master equation 
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is written in closed form for M. Thus the master equa- 
tion (5) can be expressed as a function of M. Let P(M) 
be the probability that the system has the total mag- 
netization M. The master equation for P{M) is given 

by 



dP(M, t) _ 1_ 
dt t 

N + M exp[-/?( J(M -1)/N + H)} 
2 2 cosh[/3( J(M -1)/N + H)] 
N-M exp[/3( J(M + 1)/N + H)} 



P(M) 



P{M) 



+ 



+ - 



2 2 cosh[/3( J(M + l)/N + H)] 
N - M + 2 exp[/3( J(M - l)/iV + if)] 

2 2cosh[/3(J(M-l)/AT + ii")] 
AT + M + 2 exp[-/3(J(M + 1)/N + H)\ 



P(M - 2) 



2 cosh[/3( J(M + l)/N + H)} 



P(M + 2) 



(7) 



where A" is the number of spins and M takes discrete 
values -N, -N + 2, • • • N (see Appendix A). For fi- 
nite AT we can solve the simultaneous equations for 
P(-N, t), P(-N + 2, t), ■ ■ ■ , P(N, t), as well as perform 
a Monte Carlo simulation of the model. 



III. ASYMPTOTIC SIZE-DEPENDENCE OF 
THE RELAXATION TIME NEAR THE 
SPINODAL POINT 

Next, we study the scaling properties of the relaxation 
time near the spinodal point. For large N, we put m — 
M/N and regard m as a continuous variable. Expanding 
the RHS of Eq. (7) in a series in e = 1/N, we obtain the 
following Fokkcr-Planck equation: 

dP(m) d , . d 2 , . _ o. 

dT = d^3i(m)P(m) + e—g 2 (m)P(m) + 0{e 2 ), 



where 



(8) 



(3Jm 



g 1 (m) = m-t<mh[(3(Jm + H)]+e- 

cosh [p{Jm + H)\ 

92{fn) = 1 — mtanh [(3(Jm + H)] . 

(9) 

The last term of g± (m) gives the correction to the spin- 
odal point due to the finite-size effect. Hereafter, we ig- 
nore this term as it is very small. Near the spinodal point, 
we expand gi(m) and gi (m) around the spinodal point. 
We set x = m — msp, and y — [3{H — Hsp) = h — hsp, 
and then we have 



g 1 {m)~-jfj~r 1 (pjx 2 + 2xy) 
92{m) * jj, 



(10) 



where r/ = |msp| = \J\ — 1/(3 J. 



Let us consider the time evolution of x, starting from 
x = 0. The distribution of x evolves according to Eq. (8) 
with gi ~ —y/(3J and gi ~ 1/(3 J. When x approaches 
x ~ y 1 / 2 , the correction term —r\(3Jx 2 in g\ becomes 
relevant, but other correction terms in gi, such as —2rjxy, 
are still irrelevant. Therefore, in order to determine the 
relaxation time, we can use the approximation that g\ ~ 
—y/(3J — rj(3Jx 2 in the early stage of the phase change. 
Then, the Fokker-Planck equation takes the form 



d_ 

di 



dx \(3J 



(3 J dx 2 J 

Now, we introduce the scaled parameters, 

t = *\y\- 1/2 , 



and 



A = N 2 / S y. 



(11) 
(12) 

(13) 



Then, for nonzero y, the Fokker-Planck equation is given 

by 



d_ 

di 



= < -l/3 



{-IAI 1 / 2 



+ 



d_ 

d 2 



pj\a\ oe 



P(U), 



(14) 



where A > for the upper sign, and A < for the lower 
sign. Because Eq. (14) depends only on A except for the 
factor e 1 / 3 , the relaxation time is expected to be given 
in the form 



iV 1 /3 /(A )^iVi/3 /(i v2/ 3 ^_^ p) ^ 



(15) 



This is the finite-size scaling for the relaxation time near 
the spinodal point. This form has been pointed out by 
Paul, et al.[6]. Here, we reexamine the form of the scaling 
function by studying the asymptotic forms of the relax- 
ation time. 

In the above argument leading to the scaling relation, 
it is necessary to pay attention to the range of the pa- 
rameters, in which the application of the above estimate 
can be verified. We may regard the relaxation time as a 
time when the magnetization x becomes 0(1). Then, £ 
becomes 0(|y| -1 / 2 ). This implies that we may need to in- 
clude an additional y-dependence in the relaxation time. 
Thus, we cannot immediately conclude that the finite- 
size scaling has the simple form of Eq. (15). Indeed, such 
consideration is essential for the relaxation after rapid 
quenching at zero field and has been studied as a scaling 
theory of the relaxation at unstable point [10]. In fact, 
the relaxation time is proportional to lnA^ in the relax- 
ation at the unstable point, which is unexpected from 
the form of the Fokkcr-Planck equation. This problem 
is investigated in Appendix B, where we show that the 
additional contribution does not need to be taken into 
account in the relaxation at the spinodal point, and the 
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finite-size scaling is correctly given by Eq. (15). In the 
following, we investigate asymptotic forms of the relax- 
ation time for the cases y < 0, y = 0, and y > 0. 

In order to consider more general cases later (Sec. V), 
we rewrite the Fokker-Planck equation with general co- 
efficients a, 7, and 5 : 



|^*)= £ i/3 {-iAi i/2 |(±«+7a 

5 d 2 1 

The coefficients a, 7, and 5 are defined as follows: 

gi ~ -ay - jx 2 
92 ^ S. 



(16) 



(17) 



In the case of the relaxation from the mean-field spinodal 
point, a = S = l/(/?J), 7 = r\(5J. 



Metastable region: the y < case 




FIG. 1: The scaled potential £/(£) in Eq. (19). A: the 
metastable state. B: a point in the basin of attraction of 
the stable state. C: the unstable state. 



This case corresponds to the relaxation from the 
metastable state to the equilibrium state. We can esti- 
mate the relaxation time by Kramers' method of escape 
over a potential barrier [11]. We rewrite Eq. (16) in the 
following form (A < 0) 



at 

where 



dP_ £l/3 d_ 



dZ V|A| 



tf(0=a£-7£ 3 /3 



(18) 



(19) 



is the scaled potential. We consider the escape of proba- 
bility from the valley to the outside (A — ► B) as depicted 
in Fig. 1 [11]. There, the probability current er is given 
by 



A/zJLp-\M 3/2 u(t)/s®- P \M 3/2 u(0/s 



|A| e at 



p. 



(20) 



We consider the stationary current (a =const.) and in- 
tegrate it between the two points A and B. 



,1/3 



S [, 



|A| 3/2 ^(«)/<5p]B 



|A| J* ePJW/'uw/Sf^' 



(21) 



If the number of spins N is very large, we can use the 
steepest-descent method, and we have 



,i|A| 3 / 2 (7(C) + i/3J|A| 3 / 2 C/"(C)(?-C) 2 ^ 
|A| 3/2 [/(C)/<5 



-271-5 



I A| 3/2 [/"(C) 



(22) 



This approximation is valid for a sufficiently large |A|. 
We consider the early stage of the relaxation and assume 
the relaxation has not occurred yet, namely P(B) w 0. 
Then, we obtain the estimate: 



a = P(A)e 



1/3 Z^El^lp-\ A \ 3/2 (U(C)~U(A))/S 

I 2tt|A| 1 /2 



(23) 



The probability distribution near the point A is approx- 
imately given by 



exp 



|A| 



P(0 * n A - 
where Z is a partition function 

|A|3/2 



{U{A) + \U"{A)(Z-A) 2 } 



(24) 



Z = I d£ exp 

OC 

\A\ 3 ' 2 U(A)/8 



s U{A) + -U"{A)^-Af 



2nS 



\A\V 2 U"(A) ' 



(25) 



Namely, the probability distribution near the point A is 
given by the equilibrium distribution of the approximate 
harmonic potential. The variable tia represents the total 
probability near the point A. This quantity evolves as 



1 



dt 



n A 



-n A , 



(26) 



where r is the relaxation time that we want to know. As 
the probability at the point A is given by 



P(A) ~ n A - 



exp 



I A I 3/2 

^-U{A) 



UA 



\A\ 3 / 2 U"(A) 
2~^5 ' 



(27) 



(see Eq.(24)), combining Eqs.(23), (26), and (27), we ob- 
tain 



Ar 1 /3 27r |[/"(A)Z7"(C)|- 1 / 2 |A|- 1 / 2 
■|A|3/2 



x exp 



-(U(C)-U(A)) 



(28) 



This is the result of the well-known Kramers' formula for 
the escape rate, and it agrees with the finite-size scaling 
equation (15). 
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The potential U(£) is [/(£) =a£- 7£ 3 /3, and the two 
points A and C are given by the condition dU/d£_ = 0. 
Therefore, 



A=-C=- 



a 
7 



U(C)-U(A) = ^A 



(29) 



U"{C) = -2y/ary. 
Hence, the relaxation time for sufficiently large |A| is 



ATV3 E exp /l^ 



|A| 



.3/2 



(30) 



In the case of Eq. (10), a = 5 = l/(/3J), 7 = ?7/3J, and 
we have 



ATl/3 



2 v /^ eXP l3W/ 2 



(31) 



Here, it should be noted that from Eq. (2) the following 
relation holds: 

A 3 / 2 = pAF = (3(F(C) - F(A)). 



3/?Jr? 1/2 

Therefore, the relaxation time obtained above is roughly 
r ~ e@ AF , which is the well-known Arrhenius formula. 
Another derivation of Eq. (30) uses the WKB approxi- 
mation as discussed by Tomita, et al. [12]. We can derive 
the same result (see Appendix C). 

In the infinite long-range model, microscopic fluctua- 
tions do not grow to become macroscopic because the 
long-range interaction suppresses clustering and nuclc- 
ation. It should be noted that, in the limit of N — > 00, 
the system remains at the metastablc or marginally sta- 
ble point, and the relaxation time from that point be- 
comes infinite. In contrast, in systems with short-range 
interactions, the nucleation process causes the system to 
relax to the equilibrium state in a finite relaxation time. 
Thus, the divergence of the relaxation time does not take 
place, and so far the divergence of the relaxation time 
has not been considered seriously. However, as has been 
pointed out [3], effective long-range interactions appear 
in systems in which elastic deformation mediates interac- 
tions among the spins. In such systems, the long-range 
interaction model is effectively realized, and we expect 
that the finite-size scaling discussed here would be rele- 
vant. 



At the spinodal point: the y — case 

Next, we consider the relaxation just at the spinodal 
point, y = 0. Substituting y — in the Fokker-Planck 
equation (11), we obtain 



d_ 

di 



(32) 



Putting x — e 1 / 3 z, 



dz 



By using the scaled variable s — te 1 / 3 , we can eliminate 
the e-dependence. Thus, as pointed out by Kubo et al. 
[9] , the relaxation time behaves as 



T OC £ 



-V3 = TV 1 / 3 . 



(34) 



The relaxation time diverges in the limit of N — > 00 
just at the spinodal point. In the limit of N — > 00, the 
system remains at the unstable point, and the relaxation 
time becomes infinite. This divergence is again due to 
the long-range interaction. 



Unstable region: the y > case 

Finally, we consider the case y > 0. In this case, even 
if N = 00 (e = 0), the relaxation takes place. Namely, 
the relaxation time saturates at a finite value at large 
N . Therefore, we consider only the limit N —> 00. The 
Fokker-Planck equation (14) then becomes 



dP 
~dt 



(35) 



Therefore, the relaxation time is expected to scale as 

t-jT 1 / 2 -(h-hsr)- 1 / 2 . (36) 

In the limit of N — > 00, there is no diffusion term, so we 
can derive the time evolution of the scaled magnetization 
directly. Namely, putting P(£,t) = 5(£ - £(£)), we 
obtain 



The solution of Eq.(37) is given by 



tan 



(37) 



(38) 



for yjoTfyt < tt/2. At a time \J cryyt — tt/2, the above ex- 
pression indicates that £(t) would diverge. However, the 
higher order terms in the original Fokker-Planck equation 
(8) prevent this divergence of the magnetization. Thus, 
the relaxation time is estimated as 



In the scaling form, 



(39) 



(40) 



This result is consistent with that obtained by Binder 
[1], who showed that for h > h$p, the relaxation time 
behaves as r ~ (h — /isp)~ 1//2 - 
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IV. NUMERICAL RESULTS 



of the relaxation time with the formula, 



We have seen that the relaxation time r obeys the 
finite-size scaling form (15): 

r(y,N)=N 1 / 3 f(N 2 / 3 y), 

and we have derived the asymptotic forms of the relax- 
ation time, i.e., in the case of Eq. (10), 



/(A) 



TrfalAI)- 



-1/2 



exp 



-(?7Ar 1/2 for A > 1 



(3W7i' A ' 3/2 | fOT - A>>1 

(41) 



ln(N" 1/3 x) 




FIG. 2: (Color online) Comparison of the asymptotic form of 
the relaxation time with numerical results. Data are plotted 
in the coordinates (A, ln(A r ~ 1/ ^ 3 T)). The solid lines denote 
the asymptotic form (41). The symbols are interpreted as 
follows: Closed circles, upward triangles, downward triangles 
and closed squares denote the data obtained by the master 
equation for N = 1000, 2000, 4000 and 10000, respectively. 
Open circles and squares denote the data obtained by the 
Monte carlo method for iV = 1000 and 10000, respectively. 



In this section, we check these results by numerical 
studies. We solved the original master equation (7) nu- 
merically, and also performed kinetic MC simulations 
(see Appendix D). The parameters are set as j3 = 1 and 
J = 2. The relaxation time is defined as the time at 
which the magnetization of a sample reaches a certain 
value mo. Here we adopt too = 0. In the Monte Carlo 
simulations, the relaxation time is measured directly in 
each simulation. On the other hand, in the master equa- 
tion we have to define it from the change of the probabil- 
ity distribution P(M,t). Namely, we obtain the average 



/>oo 

J° M<0 



and its standard deviation as 



/>oo 



(42) 



(43) 



M<0 



We plot data for various parameters in a scaling plot 
in Fig. 2, namely, in the coordinates (A = N 2 / 3 y, 
ln(AT -1 / 3 T)). We confirmed that both methods give the 
same results, as they should. In the Monte Carlo simu- 
lations, each data point is an average over 1000 samples, 
and the error bars are smaller than the symbol size in the 
figure. All the data collapse well onto a scaling function, 
which indicates that the finite-size scaling works well. For 
large A, data points for different N deviate slightly from 
the scaling function. This fact is explained as follows: 
The condition for the finite-size scaling to hold is that 
the system size N is sufficiently large and the magnetic 
field is sufficiently close to the spinodal point. This im- 
plies N 3> 1 and \y\ <C 1. However, an even stronger 
condition is required for the finite-size scaling. As we as- 
sumed x — to — TOgp <C 1 to derive the Fokker Planck 
equation and x was rescaled as x = (.ly] 1 ^ 2 , not only 
\y\ <C 1 but also \y\ x ^ 2 <C 1 was necessary. Therefore, 

| A |i/2 <A ri/3 
is required for the finite-size scaling. 




FIG. 3: (Color online) Ratio of the standard deviation of the 
lifetime, ay, to the lifetime, r, as obtained from the master 
equation and MC simulations. The data are shown as a scal- 
ing plot vs the scaling variable A. The symbols have the same 
interpretations as in Fig. 2. 



Now, we compare the numerical results of the master 
equation with the asymptotic form for the relaxation time 
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t, Eq. (41). In Fig. 2, we compare numerical results 
and the asymptotic forms, Eq. (41). Here we find that 
the asymptotic formulae describe the scaling form well in 
the large-A region, where the data points approach the 
asymptotic forms when the size increases. Here, we find 
that the scaling property (15) holds, and the asymptotic 
forms also hold asymptotically. 

As the scaling variable A goes from large positive to 
large negative values, the lifetime distribution goes from 
a narrow peak in the unstable region to an exponential 
distribution deep in the metastable region. This is illus- 
trated in Fig. 3 by the ratio of the standard deviation 
of the lifetime, a T , to the lifetime, r, as obtained from 
the Master equation and also by MC simulations. These 
results also show good scaling collapse, although the con- 
vergence to the asymptotic limit is somewhat slower for 
<7 r than for r, as one would expect for a quantity involv- 
ing a second moment of the lifetime distribution. 



V. THRESHOLD PHENOMENA UNDER 
EXTERNAL PUMPING 

In this section we study switching in spin-crossover 
(SC) materials under photo-irradiation, which is repre- 
sentative of a wide range of switching processes under 
external driving forces. The pumping effect is described 
by the following process in the master equation. The 
pumping operation causes a down spin (LS) to flip to 
the up state (HS) with a transition rate a per unit time. 
Thus, the total transition rate (6) from down to up is 
now given by 



w(- -> +) 



1 



-PE+ 



e -0E + + e -0E. 



+ a, 



(44) 



where E± denote the energies with the spin + and — , 
respectively. On the other hand, the transition rate from 
up to down remains unchanged and is given by 



u;(+-»-) 



1 



-/3E- 



t e-f 3E + + e-f iE - 



(45) 



These transition rates do not satisfy detailed balance, 
but they produce a stationary state. Thus, as far as we 
consider the system at a given temperature, we have a 
well-defined relaxation process to the stationary state. 
Here, we study the parameter dependence of the relax- 
ation time. To express the pumping, we add the term 

N-M N-(M-2) 
-a^—P(M) + a L Lp(M - 2) (46) 

in the master equation. In the large N limit, this modi- 
fication gives the additional terms in the Fokker-Planck 
equation 

d d 2 

~ a d^ {(1 ~ m)p(m ' t)} + ae dm^ {(1 ~ m)p(m ' t)} ■ 

(47) 



Thus gi(m) and g^m) in Eq. (9) arc now replaced by 
gi(m) = m — tanh[/3( Jm + H)] — o(l — to) 



g2(m) = 1 — mtanh [(3(Jm + H)] + a(l — to). 



(48) 



The point corresponding to the spinodal point is given 
by 



d 

gi(rn) = and —g^m) = 0, 



(49) 



which are given by 



1 + a c H 
tanhz= z - —(1 + a c ) - a c , (50) 



and 



cosh 2 , 



J3J_ 

1 + a c ' 



(51) 



where z = P(Jm sp + H). For the SC material, H consists 
of the following two ingredients: the crystal field which 
gives the energy difference between HS and LS states, 
and a time-dependent field which represents the different 
degeneracies of the HS and LS states [13]. 

Expanding the variable near the spinodal point: x = 
to — TO sp and y = a — a c , we have 



gi ~ -ay - 



where 



and 



and 



a = 



1 + a c 



l-Wi- 



1 + Q C 

(3J 



7 = /?J(l + a c )Wl- 



1 + a c 



(52) 



(53) 



(54) 



1 



92 



PJ l + a, 



' f l + Jl- 1 -^\^5. (55) 



I5J 



As gi and gi have the same forms as Eq. (17), the same 
analysis as in the previous section remains valid, and the 
finite-size scaling is again given in the form 

T = N 1 ^f{N 2 ^{a-a c )). (56) 

The asymptotic form is also given by Eqs. (30) and (40). 

VI. SUMMARY AND DISCUSSION 

Mean-field type critical behavior takes place in sys- 
tems with effective long-range interactions, as has been 
pointed out for spin-crossover type materials [3] . We ex- 
pect that the dynamical critical properties such as the 
spinodal phenomena are realized in those systems. In 
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models with short-range interactions, there exists a mode 
of relaxation from the metastable state through nucle- 
ation of localized clusters. Thus, the relaxation time 
around the spinodal point changes smoothly, and the 
critical properties at the spinodal point in the mean-field 
theory are smeared out. However, in long-range inter- 
action models, the relaxation time diverges as described 
by the mean- field theory. It is, therefore, necessary to 
study the finite-size scaling properties of the critical be- 
havior. Thus, we here studied the size dependence of the 
relaxation time near the spinodal point in the Husimi- 
Temperley model. Starting from the master equation for 
the probability density of the total magnetization un- 
der the Glauber dynamics, the Fokker-Planck equation 
for large TV was obtained as in previous work [6]. Using 
this Fokker-Planck equation, we investigated the relax- 
ation processes near the spinodal point. As a result, we 
obtained asymptotic forms and a finite-size scaling func- 
tion for the relaxation time, which cover both sides of the 
spinodal point, i.e., the metastable side and the unstable 
side. 

We further extended the analysis to systems that are 
pumped by an external field. The critical properties ob- 
tained in the present work should apply widely to thresh- 
old phenomena in long-range interacting models, such as 
the threshold phenomena found in the excitation process 
by photo-irradiation from the low-temperature phase to 
a photo-excited high-temperature phase in spin-crossover 
materials [5]. 

In the ongoing development of nano-size fabrication, 
switching processes in nano-size systems are an impor- 
tant issue, and the size-dependence of relaxation times 
of the switching processes must be precisely understood. 
We hope the scaling properties presented here will help 
to analyze such processes in experimental systems. 
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APPENDIX A: THE DERIVATION OF THE 
MASTER EQUATION 

In this Appendix, we derive the master equation (7) 
for the Hamiltonian (1) and the transition rate (6). The 
probability of the state {cti,<T2,--- ,o~n} at a time t, 



which is denoted by P(ai, • • • , (JN\t), evolves according 
to 



d 



N 



— P{<Tl, ■ ■ ■ ,a N ;t) = -^Wlf(l7i -> -<7j)P((7l, ■ ■ ■ ,GN]t) 



i=l 



N 



+ u M -2<ji (-o-j -> c7i)P(- ■ ■ , -a h •••;£) (Al) 



where 



%(<r. — ► — 0i) = 



1 exp (-(3ai(J(M - aJ/N + H)) 



r 2 cosh (fi(J(M - a t )/N + H)) ' 

(A2) 

and M is the total magnetization, i.e., M = J2i a i- We 
consider the time evolution of the probability of M: 

(71 ,<T 2 ,"- ,CT]V=±1 \»=1 / 

where 5(a, b) denotes the Kronccker delta. After some 
calculation from Eq. (Al), the equation of motion for 
P(M, t) is obtained in the form 

d N + M 

j^P(M, t)=- _^_ Wm (+1 - -1)P(M, t) 



N-M 



wjif(-l-*+l)P(M,t) 



+ - ^ 2 W 2 (-f - +l)P(M - 2,t) 

+ N+i f + 2) UM+2(+l - -l)P(M + 2,t). 

(A3) 

The meaning of this equation is clear. The first and sec- 
ond terms correspond to the transition from the magne- 
tization MtoM-2 and M + 2 respectively. The third 
and fourth terms represent the transitions from M — 2 
to M and from M + 2 to M, respectively. This equation 
gives Eq. (7). 



APPENDIX B: CUT-OFF INDEPENDENCE OF 
THE FOKKER-PLANCK EQUATION (16) 

In Sec. II, we remarked on the possibility of an ad- 
ditional y-dependence in the relaxation time. Namely, 
if we regard the relaxation time as the time when the 
magnetization x becomes O(l), this corresponds to the 
time when £ becomes 0(|y|~ 1//2 ), and this implies that we 
cannot conclude the finite-size scaling of the relaxation 
time, Eq. (15), from the form of the Fokker-Planck equa- 
tion (16). In other words, although the Fokker-Planck 
equation (16) seems to depend only on A, we must re- 
strict the range of the variable £ < 0(y~ 1//2 ) and this 
cut-off of £ can induce an additional y-dependence in the 
relaxation time. We note that the relaxation time indeed 
depends on the cut-off in other situations. One exam- 
ple is the relaxation from the mean-field unstable fixed 
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point. In this case, the Fokker- Planck equation is given 
by 

Wt P{x ' t) = ~^ xP{x ' t) + e ^ P{x ^ (B1) 

(we set some coefficients equal to unity). We can trans- 
form this equation to the scaling form similarly. If we set 



: 1/2 C, 



(B2) 



This equation is apparently independent of e. Is the re- 
laxation time independent of the system size N = 1/el 
The answer is No. It is known that the relaxation time 
in this case is t ~ hxN [10]. We show that this TV- 
dependence stems from the finite cut-off. We can solve 
Eq. (B2) for the initial condition P(f,0) = <5(£), 



where a(t) is given by 



2na(t) 



exp 



e 



2a{tf 



cr(t) = y/# 



1 



(B3) 



(B4) 



It takes infinite time for £ to reach infinity. As x = e 1 ' 2 £, 
x 2 (t) ~ ea(t) 2 = ee 2t . We consider the relaxation time 
as the time when x 2 (t) reaches 1, i.e. x 2 (t) ~ 1, the 
relaxation time is proportional to In TV, 



In N. 



(B5) 



In this way, we found out that the cut-off dependence 
could actually affect the relaxation time, but this cut- 
off played no role in the case of the relaxation near the 
spinodal point. 

Hence, here we show that the relaxation time does not 
depend on the cut-off of £ if this cut-off is very large. 

If we denote the average of £ n over P(£)by £™, the time 
evolution of £ is given by 



=p 1 /3| 



^e^lAI 1 / 2 (±a + 1 at) 2 ) 



(B6) 



If £(io) is larger than yj a/7 = v, it can be shown from 
Eq. (B6) that 



£(*)> 



where A(t) is given by 



A(t) = 



e(*o) 



■ exp 



1 + A(t) 
l-A(t)' 



2s 1 / 3 h" 2 vl (t - 1 ) 



(B7) 



(B8) 



The average of the scaled magnetization £ (t) reaches in- 
finity when A(t) = 1, namely 



t-tn 



hi 



e(*o) 



(B9) 



Because to is finite, it takes only a finite time for £(t) 
to reach infinity. Therefore, there is no cut-off depen- 
dence on the relaxation time in the Fokker-Planck equa- 
tion (16). 



APPENDIX C: DERIVATION OF EQ. (30) BY 
THE WKB APPROXIMATION 




A 

metastatic 



equilibrium 




equilibrium 



FIG. 4: Rough sketches of the scaled potential U(£) and the 
corresponding Schrodinger potential V(£). 



In the body of this paper, we estimated the relaxation 
time for A < and |A| 3> 1 according to Kramers' ar- 
gument. Here we give another derivation by using the 
WKB approximation. The following derivation is es- 
sentially the same as that of Tomita, et al. [12]. The 
Fokker-Planck equation (18) can be transformed to the 
"Schrodinger equation" 



dQ = 1/3 S d 2 Q 

ot ~ e |A| &e 

by substituting 



: 1/3 V(0Q(0=e 1/3 HQ(0 (CI) 



r = «'x,> ( -1|A| 3 / 2 e/(o) Q(0- 



The scaled potential £/(£) is given by Eq. (19), and the 
Schrodinger potential V(£) is 

V(0 = ±\A\ 2 U'(0 2 -l\M 1/2 U"(0- (C2) 

If the eigenvalues of 7i are \ (i = 0, 1, 2, • • • ) and the 
eigenf unctions are <pi, we can expand Q(£,t) as 



r 1/3 A,t 



(C3) 



The lowest eigenvalue is Ao = 0, and the corresponding 
eigenfunction is 



<M0 



1 



z 1 / 2 



-\A\ 3/2 U(t)/(25) 



(C4) 
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which corresponds to the equilibrium state. Z is the 
normalization factor for j </>g<i£ = 1. The second lowest 
eigcnfunction <fii will represent the metastable mode, and 
the corresponding eigenvalue will be connected with the 
inverse of the lifetime of the metastable state, e 1 / 3 Ai ~ 
1/r. 

From the Schrodinger equation (CI), 



Substituting Eqs. (C12) and (C13) into Eq. (C8), 



Q = - 7 L<pi> + V ( Po 
2m 

Ai^i = --z—<P'{ + V(j>i. 
2m 



(C5) 
(C6) 



The mass m is m — |A|/(2<5). If we multiply Eq. (C5) by 
4>\ and Eq. (C6) by tfio and subtract the two equations, 
we obtain 

Ai^i = ^(^i-Mo)'- (C7) 

Integrating this equation from — oo to the point C (see 
Fig. 4), we obtain 



Ai 



2m I-oc <h<t>id4 



(C8) 



because <Pq{C) = 0. Here we consider the metastable 
wave function <fi ms , which corresponds to the localized 
canonical distribution at the valley of the potential. 
Hence we assume for £ < C 



0o ~ u<f) ms for £ < C, 
where u is a constant given by 



J —oo 



2tt5 



(C9) 



\A\ 3/2 U(A)/S 



\A\ 3 / 2 U"(A) 

(CIO) 

Besides we assume that the first excited eigenfunction is 
also proportional to </> ms in the range £ < C, 



h ~ w^ms for £ < C, 



(Cll) 



because 4>i is considered to represent the metastable 
mode. Under these assumptions, we get 

/C nC 
0o<M£ = uv 4> 2 ms dt » uv. (C12) 
-oo J —oo 

Using the WKB approximation, it is obtained that 



b[(C) w ^2mV(C)MC) - -V2mV(C)MC). (C13) 



Ai — 



VjC) 
2m 



(C14) 



After some calculation, wc obtain 



*i = U—\U"(A)U"(C)\e-P AF (C15) 

2 V 7T 

where the free energy barrier [3AF is /3A_F = 
|A| 3 / 2 (J7(C) - U(A))/5. Therefore the lifetime of the 
metastable state r, which is equivalent with the relax- 
ation time, is 

r ~ 2N 1 /^\U"{A)U"{C)\- 1 ' 2 \K\- 1 ' 2 el' AF . (C16) 

Comparing with the relaxation time obtained by 
Kramers' argument, Eq. (28), they agree with each other 
except for the minor difference in the constant prefactor. 



APPENDIX D: MONTE CARLO SIMULATION 

We also obtained data by performing kinetic Monte 
Carlo simulations to confirm the data obtained by solving 
the master equation. In the Monte Carlo simulations, 
each data point is an average over 1000 samples, and the 
error bars are smaller than the symbol size in Fig. 2. 

The algorithm of the Monte Carlo simulations is as 
follows. We choose a spin at a site i randomly, and up- 
date the spin with the probability corresponding to the 
Glauber model given by Eq.(6): 



_ cxp (-P<Tj(J(M - <7 % )/N + H)) 
Gt) 2 cosh (P(J(M -<Ji)/N + H)) ' 



where ai is the spin state of the i-th spin (<t, = ±1). In 
principle, a small time increment At <C 1 is necessary 
to reproduce the result of the master equation given as 
a differential equation. However, we found almost the 
same result with the different time division, A£=0.01 and 
1 for the quantities plotted in Fig. 2. Hence we obtained 
the data with At = 1. During one Monte Carlo step we 
perform single spin flips N times. Therefore, the time t is 
related to the Monte Carlo step s by t = s because At = 
1. The initial condition of each Monte Carlo simulation 
is set to the spinodal magnetization. 
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